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DEVELOPMENT AND CORRELATION: VIKING ORBIT ER 
ANALYTICAL DYNAMIC MODEL WITH MODAL TEST* 


B. K. Wada, J. A. Garba, and J. C. Chen 
Jet Propulsion Laboratory 
Pasadena, California 


The Jet Propulsion Laboratory is responsible for the Viking Orbiter 
System, which is part of the overall Viking Project managed by the 
Viking Project Office at Langley Research Center for NASA. 

The development of a mathematical dynamic model and its verification 
by a modal test is a significant milestone for many Projects including 
Viking Orbiter (VO). Difficulties encountered include performing a 
modal test, establishing a criteria for correlation of analysis with test, 1 
and modifying a large finite element mathematical model to match test 
data if required. Often the modal test is performed near the end of the 
Project development schedule; consequently, the time alloted to obtain 
a verified mathematical model is minimal. 

The paper describes the VO experience in the achievement of a good 
mathematical model. Success can be attributed to the coordination of 
analysis and tests using substructure modal coupling techniques. The 
experience would benefit the overall planning of any project, such as 
Shuttle, especially if substructure modal coupling techniques are 
contemplated. 


INTRODUCTION 

The Jet Propulsion Laboratory (JPL) is 
responsible for the Viking Orbiter System, 
which is part of the overall Viking Project 
managed by the Viking Project Office at 
Langley Research Center (LRC) for NASA. The 
Spacecraft will be launched on a Titan IIIE/ 
Centaur Launch Vehicle in August 1975. 

The total launch vehicle system consists 
of numerous subsystems that are developed by 
various aerospace organizations. The creation 
of a launch vehicle system model requires the 
transfer of each organization's mathematical 
models. One organization eventually creates 
the total model required for analyses. The 
complexity and size of the problem required 
the use of substructure modal coupling concepts. 


To minimize schedule and cost, the goal was to 
limit the responsibilities of each organization 
to their own mathematical models and verifi- 
cation test program. A strong emphasis on 
technical accuracy existed. 

The development of a test-verified mathe- 
matical dynamic model is a significant mile- 
stone for many projects including Viking 
Orbiter (VO). Difficulties encountered include 
performing a modal test, establishing a crite- 
rion for correlation of analyses with the test, 
and modifying a large finite element mathemat- 
ical model to match test data if required. 

Since the modal test was performed near the 
end of the Project development schedule, the 
time available to obtain a test -verified mathe- 
matical dynamic modal was minima.1. This 
paper describes the VO plans and experience to 
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obtain a good model. Emphasis was placed on 
the early development of good mathematical 
models, performance of the modal test, and 
methods to correlate the analysis with test 
data. 

A valid mathematical model for VO was 
required because the design and flight loads 
for the primary structure were established by 
load analysis. Load analysis is a procedure 
for obtaining VO member forces from the 
dynamic response of a complex finite element 
model of the complete Launch Vehicle System 
{including the VO) subjected to launch vehicle 
engine transients. Load analysis requires sub- 
structure modal coupling (Ref. 1) of the various 
structural subsystems of the Launch Vehicle 
System to allow a solution within present com- 
puter capabilities. In addition to the usual 
objectives of modal tests, determination of 
individual member forces is emphasized 
throughout the program. 

A good mathematical model was generated 
by establishing an overall plan integrating sub- 
system analysis and test with the substructure 
modal coupling approach. Thus the model was 
continually updated during the program. 

Emphasis was also placed on the modal 
test and the establishment of a measure of 
correlation of the analysis with the test. The 
measure of correlation is required to establish 
a factor directly related to the confidence 
placed in the member forces resulting from 
load analysis. 

This paper describes three general activi- 
ties that resulted in the VO analytical dynamic 
model, and that were updated and verified by 
test data during the Project. 

(1) The generation of the overall plan for 
load analysis, an analytical dynamic 
model, and development tests. 

(2) The performance of VO subsystem 
static and modal tests. 

(3) The correlation of the VO System 
modal analysis and test. 

The details of the modal test are not 
included (Ref. 2). However, actual results are 
used to show the degree of success attained on 
a large complex structure. The substructure, 
tests, and update of substructure mathematical 
models occurred between July 1, 1972 and 
May 15, 1973, and the VO System modal test 
between June 1, 1973 and July 30, 1973. The 


correlation of modal test results with analysis 
occurred between July 1, 1973 and July 30, 
1973. The final mathematical model was 
completed on schedule by July 30, 1973. 


DESCRIPTION OF HARDWARE 

Figure 1 identifies parts of the Viking 
spacecraft (V-S/C), Viking transition adapter 
(VTA), and Centaur truss adapter (CTA) per- 
tinent to this discussion. 



Fig. 1 - Viking spacecraft 

The Viking Orbiter System was compli- 
cated because it is situated between the Viking’ 
lander capsule (VLC) on top and the Centaur 
adapter (VTA/ CTA) on the bottom. The VLC 
and the VTA/C TA are the responsibility of 
Martin Marietta Aerospace (MMA) and General 
Dynamic s/Convair Astronautics (GD./CA) 
respectively. 

The weight of the hardware is summarized 
in Table 1. 

Figure 2 is a description of the VO/VTA/ 
CTA configuration for the modal test per- 
formed at JPL. 

Figure 3 is a photograph of the test of the 
orbiter development test model (ODTM). The 
rationale of the configuration will be discussed. 
The major differences between the test and 
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TABLE 1 

Approximate Weights* 


Hardware 

Weight 

(lb) 

Responsible 

Organization 

VLC 

2567 

MMA 

VOS 


JPL 

Bus and adapters 

1109 


High -gain antenna with support 

47 


Scan platform with support 

201 


Four solar panels 

244 


Cable trough 

49 


Propulsion module hardware 

512 


Propellants*! 

3138 


VTA 

t 

GD/CA 

CTA 

X 

GD/CA 

Total weight 

7867 

LRC 

*The weights are the values used for analysis on 7/1/73. 
tSee Table 2. 

^Included in the Centaur Model. 


TABLE 2 

Approximate Propellant Weight and Ullage Summary 



Oxidizer 

Fuel 

Configuration 

Fluid 

Rigid 

Weight 

(lb) 

Slosh 

Weight 

(lb) 

Ullage 

(%) 

Fluid 

Rigid 

Weight 

(lb) 

Slosh 

Weight 

(lb) 

Ullage 

(%) 

Viking Mission A1 

N 2°4 

1379 

431 

20.9 

MMH* 

968 

235 

13.3 

Viking Mission A2 

N 2°4 

1318 

428 

23.7 

MMH 

911 

249 

16.36 

Viking Mission B 

N 2°4 

1470 

415 

17.62 

MMH 

1049 

204 

9.68 

Modal test 

Freon- 

TF 

1735 

371 

17.62 

Alcohol 

902 

227 

9.68 

*Mo no methyl hydrazine. 
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NUMBERS IN PARENTHESES ARE NODE NUMBERS 



Fig. 2 - Test configuration and 
node identification 


flight configurations are that the test configura- 
tion has: 

(1) Rigid VLC with inertia property 
simulation. 

(2) No solar panels and solar panel 
dampers. 

(3) No high gain antenna. 

(4) Propulsion propellant mass loading as 
shown in Table 2. 

(5) Mass and stiffness to represent 
dynamic characteristics up to 60 Hz. 

(6) No slippage of the scan platform joint 
along the serrations. 

The propellant loading for different con- 
figurations is summarized in Table 2. The 
information will be of value for future 
discussion. 


GENERAL APPROACH 

The analysis plan and modal test approach 
were closely integrated with VO Project plans 


and requirements (Ref. 3). Our belief is that 
a successful development of a mathematical 
model correlated by test is directly related to 
the overall analysis, hardware, and test plan. 

A. Load Analysis and Its Impact 

The design and flight loads for the pri- 
mary structural members were established by 
load analysis. The load analysis is a dynamic 
analysis procedure to obtain VO member 
forces. The complete Launch Vehicle System 
including the VO is subjected to launch vehicle 
engine transients measured from past flights. 

The use of load analysis for design and 
flight loads necessitated a continual reitera- 
tion of the VO mathematical model to update 
the design loads as the design evolved. A VO 
model of the final configuration verified by a 



Fig. 3 - ODTM modal test 
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system test was required". The flight loads 
were used to establish forces for which the 
structure was qualified. 

The establishment of flight loads and the 
structural qualification test program were near 
the end of the development schedule. Since the 
modal test was also near the end of the develop- 
ment schedule, confidence had to exist in the 
mathematical model when the modal test was 
performed. 

Correlation of the modal test results to 
analysis was necessary to help establish a 
measure of uncertainty, which was required to 
establish the accuracy of flight loads. The 
measure of uncertainty is defined as the load 
analysis factor (LAF). The flight loads are the 
calculated loads times the LAF. 

B. Member Loads 

The significant parameter in load analysis 
is the member forces. The goal was to obtain 
accurate dynamic member forces, not accelera- 
tions. The modal test included the measure- 
ment of modal force coefficients. 

C. Substructure Modal Coupling 

The modal coupling of the VO with the VLC 
and VTA/CTA required consideration of: 

(1) Accurate selection of displacement 
functions. 

(2) Simplification of interfaces between 
organizations. 

(3) Provision for each organization to 
perform analyses and tests independent 
of the others. 

(4) Ability to verify the model by modal 
test. 

(5) Availability of test hardware. 

The structure below the CTA was modeled 
as a planar structure (plane before deformation 
remains a plane after deformation) and the 
structure above the VTA was a three- 
dimensional model. The requirement to 
modally couple the V-S/C to the Centaur 
resulted in the decision to include the VTA,/ 

CTA with the VO. Otherwise, the number of 
compatibility relationships would have increased 
along with the possibility of erroneous data 
caused by round-off errors. Consequently, the 
CTA/VTA was included as a part of the VO 


modal test. Another goal was to simulate the 
interfaces since truss joints were of concern. 

The substructure modal coupling analysis 
techniques were also used on VO to: 

(1) Provide a cost effective solution. 

(2) Allow use of substructure test data as 
available. 

(3) Decrease the effort to update the 
mathematical model based on the test 
data. 

(4) Increase confidence in the final model. 

A substructure was defined as being com- 
patible with: 

(1) Deliverable hardware used to obtain 
test data incorporable into the models. 

(2) Ease of interface definition and 
analysis. 

(3) Area of engineering responsibility. 

Substructure tests were used to verify and 
adjust the mathematical models. Errors were 
minimized since the responsible engineer of a 
substructure used engineering judgement to 
verify his mathematical model. 

D. Rigid VLC 

Various methods of modal coupling of the 
VO to the VLC were possible. The inclusion 
of the rigid VLC is mathematically equivalent 
to mass loading (Ref. 4) the VO interface with 
the VLC. The disadvantage of including VLC 
data into the VO analysis and test were offset 
by: 

(1) The capability to modify the rigid VLC 
inertia properties after the VO model 
was delivered to MMA. 

(2) The similarity of the VO displacement 
functions to V-S/C functions, thus 
fewer modes were required. 

(3) The capability to use the resulting 
representative configuration for the 
sine vibration tests. 


ANALYSIS 

The equations are developedito briefly 
illustrate the methodology in the creation of the 
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VO mathematical model. Also they are used as 
a basis to definitize the objectives of the tests 
and to define the data used to correlate the 
analyses with test data, 

A. Substructures 

The two general equations for the sub- 
structure are: 



[m]Ju‘|+ [c])u(+ [k] |u| = jOj (2) 
|P| = [S]M (3) 

where 

DO = stiffness matrix 
[c] = damping matrix 
[m] = mass matrix 


N 

u = || u f T [k]|u| = |^|u[ j T [ k ].| u |. 

j=l 

N 

D = yW T [ c ]jMj 

j=i 

(5) 

F. = force 
] 

N * number of finite 
elements 

{uLjuf [m].,[k] = parameters 

J J J J associated with the 
jth finite element. 

Other parameters for correlation of the 
analysis and test data are the kinetic energy T 
and the potential energy U. The dissipation 
function D cannot be used for correlation 
since the test data are used in the analysis. 
Equation (1) can be written as: 

[ u o] = [ k oo] ([ f o] ■ [ k oi] [ u l] ) (6) 


[f] = force matrix 
|P| = member forces 

[s] - force coefficient matrix 

(u} = displacement 

I = subscript representing interface 
degree of freedom 

O = subscript representing degree of 
freedom other than the interface 

Equation (2) can be derived from the Lagrangian 
equation: 


d /8L\ 9L_ dD 

dt j Du. + j 

\ au.f 1 Su. J 

' V i J 


L = T - U 

N 


and Eq. (2) as: 

[m]|u}+ [k]Ju| = |0} (7) 

where experimental modal damping values are 
used. The [c*] is assumed to be of a form 
where the transformation formed by the eigen- 
vectors of Eq. (7) uncouples Eq. (2). 

B. Displacement Functions of Substructures 

Often the dynamic characteristics of sub- 
structures are represented by a finite number 
of displacement functions to reduce the number 
of independent variables. The various forms of 
identification of displacement functions will be 
discussed. 

C. Rigid Body Modes 

Rigid body modes represent the motions 
[0R]°f the substructure when a degree-of- 
freeaom [uj]is displaced an arbitrary value 
without force. The is a solution to 
Eq. (6), where [fo] = 0, [uj] is a unit matrix 
in the degrees of freedom associated with the 
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rigid body modes. The displacements of the 
substructure due to rigid body modes are 

{ u r} = } 


Attachment modes are displacements 
[0 A ] of the substructure corresponding to con 
centrated loads [f A ] on the substructure. Dis 
placement [0 a] are solution to Eq. (6), 
where [ fo] = [^a! an ^ [ u l] ~ 0. The displace- 
ments due to attachment modes are 


The number of rigid body motions may range 
from 1 to =o . Rigid body motions in excess of 
6 are related to linkages within the substruc- 
ture. 

If the displacement at the substructure 
interface gridpoints can be represented by a 
linear combination of rigid body modes, the 
interface is defined as statically determinant. 

D. Constraint Modes (Ref. 1) 


Constraint modes represent the motions 
[0 C ] of the substructure when the displacement 
of an interface degree of freedom (DOF) 
requires force as the other interface degrees 
of freedom are restrained. Constraint modes 
are used to define displacement functions cor- 
responding to interface distortions. The con- 
straint modes [0c ] are the solution to Eq. (6), 
where [fg] = 0, [uj] is a unit matrix in the 
degrees of freedom associated with constraint 
modes. A force matrix [fc] associated with 
constraint modes exists. Note that a distinction 
between rigid body modes and constraint modes 
is not required. The displacements of the sub- 
structure due to constraint modes are 


{“aH*a]M (10) 


A disadvantage of attachment modes is 
nonorthogonality to the normal modes or to 
each other. Thus, unless extreme care is 
exercised, attachment modes that are nearly 
linear combinations of normal modes or other 
attachment modes may inadvertently be 
selected. If the system equations comprising 
substructure modes are not independent, the 
equations cannot be solved. 


F. Normal Modes 


The normal modes of the substructure are 
evaluated from Eq. (7). The displacements due 
to normal modes are 


r t 

{ u n1 = [v|K} 


(11) 


G. Total Displacement Function 


{ u c} M w 


The displacement of the substructure can 
be any combination of displacement functions 
selected above. 


They are defined only if the interfaces are 
statically indeterminant. (Interface cannot be 
defined as a linear combination of rigid body 
modes.) 

Features of constraint modes include 
orthogonal to normal modes evaluated with all 
interface degrees of freedom constrained. 

E. Attachment Modes (Ref. 5) 

The number of displacement functions 
necessary to represent the dynamic character- 
istics of the system may be minimized by the 
selection of substructure modes that closely 
represent system modes. In a combined struc- 
tural system, a gridpoint of a substructure at 
which another substructure is attached is sub- 
jected to concentrated attachment forces. These 
forces result in a displacement function that 
may have to be represented by many normal 
modes. 


M= [MM MM 



wM 

( 12 ) 


H. Generalized Coordinates 

Substitution of Eq. (12) into Eq. (7) and 
premultiplication by [0] T result in 

[>] T [m] [*]|q[ + M T [k] 0]|q( ■ { 0} 
or 

[m q ]{qf + [k q ]{qi= {0} 

(13) 

[S]M|q( = |P( 
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Damping corresponding to the displacement 
functions is introduced into Eq. (2) if desired. 

[m q ]{q} + [ c q ]{q} + [k q ]{q}={°} (14) 

I. Mixed Coordinates 

Equation (14) represents the equation of 
motion of any substructure in terms of its 
generalized coordinates. The VO analysis 
approach uses hybrid coordinates where the 
real displacements of the bus are retained and 
the generalized coordinates of the substructure 
attached to the bus are used. 

Mixed coordinates were used because of 
the order of the system equations of motion. 
Originally generalized coordinates for all sub- 
systems were retained as described in Ref. 6. 
This approach resulted in numerical difficulties 
due to limitations of single precision arithmetic 
in the Structural Analysis and Matrix Interpre- 
tive System (SAMIS) on the Univac 1108 com- 
puter. The use of real coordinates for all 
subsystems was rejected because of size limi- 
tations in the eigenvalue routine. The symbols 
used are shown below. 



The terms are defined. Subscripts differ- 
entiating substructures are introduced at this 
time. 

jq^ = generalized coordinates of i n 
substructure 

jufg = real displacement of the VO bus 

jujg. = subset of |u|j3 defining compati- 
bility of the ith substructure to 
the bus. Often the displacements 
are in local coordinates. 

The Eq. (14) for the i^ substructure with 
the interface degrees of freedom is 



Two classes of interfaces are distinguished. 
Statically determinant and statically indeter- 
minant interfaces result in different forms of 
the coefficient matrices [mq]p [cq]^ and 
[kq]i- 

A statically determinant interface exists 
when the interface coordinates jujgi represent 
linear combination of rigid body displacements. 
Although it can be greater than six, six is the 
maximum number for VO. The matrices are 
of the form 


RR RE n 


" 


- - 

m . m . 
qi qi 

ER EE 
m . m . 

_ qi qi _ 

9 

0 0 

n EE 
0 c . 

L v J 

, and 

0 0 
0 k EE 

V J 


where superscripts 

R = rigid body motion 
E = elastic motion 

Statically indeterminant interface coordi- 
nates {u|bi have more interface coordinates 
than can be represented by a linear combination 
of rigid body displacements. The matrices 
are of the form 


r m II 

m . 
qi 

- o 


r n 

c I0 ‘ 


fk 11 

k I0 ‘ 


qi 

qi 

, and 

qi 

qi 

m 01 

m . 

L v 

m 00 
m . 

J 

> 

c 01 

c 00 

k 01 

k°° 


fv- 

C qi J 


qi 



(17) 


where superscripts 

I = interface motion 

O = other than interface motion 

Equation (16) is a special form of Eq. (17). 

J. Equation of Motion of Total Structure 

If each substructure is considered as a 
finite element, they can be combined with the 
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where 


bus in a way similar to Eqs. (2) and (4). The 
results are 


j 

U R 

( 

\ 

Mr 

1 

< \ 

i iuiR 


| Ut B 

\ 

|Mb 

1 

He 

1 — 1 

3 

1 1 

Wi. 

+ 

1 — 1 

0 

i — 1 

Wij 

+ 

1 — t 
1 1 




1 


1 

f'J 


1pI b - C s 1! u 1b 


|p!i - [s] [^K 

where 

|u} = rigid body degrees of freedom of 

R the V-S/C or the motions at the 
V-S/C /Centaur interface 

th 

[0] = displacement function of i 

i substructure 

th 

|u[, = real displacements of i 
1 substructure 

jp} B = member forces in the bus 

th 

jp(. = member forces in 1 substructure 


The eigenvalues and eigenvectors of 
Eq. (18) with terms associated with [c] and 
ju} R removed are and [#]. Substitution 

of the transformation 


\W!i ( 

{ . W{x} < 19 > 


T 

into Eq. (18) and pre multiplication by [#] 
results in 



[Mj = 


r M RR M RE 

m er |/m ee -]. 


W=w t mw 


[c] = 


0 ["c E £ 


J 


_"c E ^]= 2 p.V M EE K 


=EE 


(21) 


where p\ is the critical viscous damping ratio, 

ro 0 


[K] = 


0 


0 E! U 


['s e O-W t Ck]W-[' sEE <| 

The significance of the rigid VLC is 
illustrated. Since the VLC attachment to the 
VO is statically determinant, from Eq. (16) 


.„RR ..RE 1 
M y M _ 
qL qL 

..ER ..EE 

M T M T 

qL qL 


0 0 
EE 


0 C 


qL 


, and 


0 0 
.EE 


0 K 


'qL 


where i = L is the symbol used for the VLC. 
Combination of the mass matrix of the VLC 
with the bus can be shown as 


r 

1 


,_EE I 

V 


M 


ER 

qL 


M 


RE 

qL 


w rr 

M qL + 


MASS OF 1 

BUS AT I 

L I INTERFACE q 

MASS OF VO 
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The damping and stiffness matrices do not 
couple. Thus the physical addition of the rigid 
VLC is equivalent to adding [niq^ ] to mass 
matrix of the VO. No other VLC parameters 
are included in the VO analysis. 

The equations are derived for a system 
with small damping where the eigenvectors of 
the undamped equation can be assumed to 
diagonalize the original damping matrix. This 
is assumed for VO although discrete viscous 
dampers exist. See Ref. (6) for treatment of 
viscous dampers. 

K. Size of the Problem 

The approximate size of the VO dynamic 
model used for Load Analysis is summarized 
in Table 3. 


SUBSTRUCTURE TESTS AND 
ANALYSIS CORRELATION 

During the program, information on sub- 
structures, structural components, and param- 
eters related to dynamics were obtained during 
the development test program. The tests were 
run to 

(1) Directly obtain dynamic data. 

(2) Establish feasibility of future tests. 


The principal objectives of many tests 
were other than obtaining data to verify models. 

Two tests to be discussed illustrate the 
procedure used in correlating the mathemati- 
cal model with the test data. 

A. Propulsion Module Modal Test 

Figure 4 shows the propulsion module 
modal test setup. The objectives of the test 
were to: 

(1) Establish the difficulty of performing 
a modal test with ullage in the tanks. 
(Both a zero ullage and ullage condi- 
tions were tested.) 

(2) Establish nonlinearity of the system 
with excitation force. 

(3) Establish the influence of the tank 
pressure on its dynamic character- 
istics. 

(4) Measure the eigenvectors and 
eigenvalues. 

(5) Indirectly measure the constraint 
modes. 

(6) Measure the modal force transforma- 
tions. 

The changes in the model as a result of the 
test are shown in Tables 4, 5, and 6. 


TABLE 3 

Size of VO Dynamic Model 


Substructure 

Elastic DOF 

Dynamic DOF 

Interface DOF 

Normal Modes 

Rigid lander and bus 

1,720 

153 

75 

0 

Scan platform 

580 

84 

14 

3 

Solar panels 

3,444 

452 

28 

20 

Cable trough 

192 

153 

20 

0 

Propulsion module 

695 

78 

16 

12 

3 -hole tab 

20,000 

0 

192 

0 

Mickey mouse tab 

3,400 

0 

24 

0 

Siamese tab 

1,760 

0 

22 

0 

CTA/VTA 

42 

0 

36 

0 
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Fig. 4 - Propulsion module modal test 


B. Solar Panel Modal Test 

Figure 5 shows the solar panel modal test 
setup. A brief description of the results is 
shown in Table 7. 

Since the mode shapes of the analysis and 
the test compared well, only the frequencies of 
the modes were changed. Extraneous modes 
from the analysis were eliminated. The modal 
test was performed on the solar panel with the 
relay antenna. The solar panel model result 
without the relay antenna was obtained analyti- 
cally using the model with the relay antenna 
adjusted to the test data. 

SYSTEM MODAL TEST AND 
ANALYSIS CORRELATION 

A. Modal Test Configuration 

As illustrated in Fig. 2, the VO system 
modal test configuration did not duplicate the 


flight configuration. The objective was to verify 
as many significant substructures and their 
interactions as feasible. The modal test con- 
figuration is represented by the equations 


[0] - [M] 


T,A 


( * U *R 


' |{,| R ) 

) 1U| B j 


|! u I b J 

\ ^P,T| 

[C] < 

I4l p T / 

P q <S,T 

T,A 

|q| S)T ( 

{^C J 


l4( c ] 




T,A 


' <R 
' U <B 
| q }p )T 

| q ( s>T 

f q t C 


Fig. 5 - Solar panel modal test 
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TABLE 4 

Frequency Change and Damping of Propulsion Subsystem Modal Test 


Mode No. 

Pretest Model 
(Hz) 

Modal Test 
(Hz) 

Posttest Model 
(Hz) 

Mode Description 

Damping Ratio 

1 

12.11 

12.95 

12.95 

Oxidizer and fuel tanks 
in Y direction (in phase) 

0.010 

2 

14.90 

17.66 

16.36 

Oxidizer tank in 
Z direction 

0.0049 

3 

19.24 

20.80 

19.83 

Oxidizer and fuel tanks 
in Z direction 

0.0084 


22.43 



Local thrust plate 


4 

25.76 

22.97 

26.70 

Oxidizer and fuel tank 
in Y direction (out of 
phase) 


5 

27.50 

28.33 

28.30 

Fuel tank in Z direction 

0.0074 


35.57 



Local PCA 


6 

38.67 

32.76 

34.25 

Pressurant tank in 0 y 
direction* 

0.0106 

7 


42.80 t 

45.75 

Local PCA in 
X direction^ 


8 

41.62 

50.67 

49.69 

Pressurant tank in 
Y direction 

0.0078 

12 

50.54 



65.38 

65.67 

PCA in direction 

0.0107 

* 

0 y is rotation in radians 
' MMA test data 


{qfg ^ = generalized coordinate of the scan 
’ platform in test configuration; 

joints are not allowed to slip along 
the serration 

jqf = generalized coordiriate of the 
cable trough 

In steps identical to obtaining Eqs. (19) through 
(21), one obtains from Eq. (22): 

eigenvalues, 

M T ,A < 22a > 

eigenvectors, 

M TiA (22b) 


E®] t u tg 

tPli ■ [®] Mi !“li 

where 

T,A = subscript representing analytical 
estimate of the modal test 
configuration 

}u|^ - rigid body displacement 

|u( B = displacement of bus 

|q{p T = generalized coordinate of pro- 
’ pulsion module with test propel- 
lant mass (see Table 2) 
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TABLE 5 

Kinetic Energy Distribution of Propulsion Module Modal Test 


Mode 

Pretest Analysis (%) 

Modal Test (%) 

Posttest Analysis (%) 


Oxidizer tank 

X 

* 





Y 

63.38 

76.76 

68.17 



Z 




12.95 Hz 

Fuel tank 

X 






Y 

21.41 

18.04 

24.27 



Z 





Oxidizer tank 

X 

33.11 

23.63 

34.02 



Y 






Z 

25.39 

50.09 

35.12 

17.66 Hz 

Fuel tank 

X 

15.65 

13.81 

16.72 



Y 






Z 

2.89 

0.71 

2.22 


Oxidizer tank 

X 

13.83 

26.96 

19.27 



Y 






Z 

58.56 

38.83 

49.82 

20.80 Hz 

Fuel tank 

X 

2.33 

6.22 

3.52 



Y 






Z 

19.03 

22.34 

22.04 


Oxidizer tank 

X 






Y 

22.14 

21.82 

23.26 



Z 




22.97 Hz 

Fuel tank 

X 

0.84 

0.02 

1.22 



Y 

50.75 

64.18 

50.08 



z 

3.27 

0.09 

3.66 


Oxidizer tank 

X 

11.60 

11.72 

11.29 



Y 

1.65 

0.37 

2.04 



Z 

5.22 

3.54 

5.75 

28.33 Hz 

Fuel tank 

X 

6.78 

8.35 

6.39 



Y 

2.74 

0.83 

3.42 



Z 

57.10 

58.92 

56.58 
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TABLE 6 


Modal Forces of Propulsion Subsystem Modal Test 


Member No. 

Mode 1 (lb) 

Mode 2 (lb) 

Mode 3 (lb) 

Mode 4 (lb) 



Pretest 


133.9 

-220.4 

-12.99 

-43.4 

4 

Test 


77.1 

-266.1 

-53.3 

-43.7 


Posttest 

-572.9 

112.9 

-237.3 

-110.5 

-41.9 


Pretest 

-70.8 

-130.3 

63.2 

135.7 

95.4 

3 

Test 

-66.8 

-87.1 

85.3 

90.5 

73.7 


Posttest 

-73.8 

-109.7 

72.7 

124.1 

97.5 


Pretest 

469.7 

-241.9 

-179.2 

231.6 

-32.0 

41 

Test 

476.6 

-252.5 

-109.6 

188.0 

-57.4 


Posttest 

427.2 

-244.3 

-147.8 

183.4 

! -30.9 


Pretest 

237.6 

158.6 

110.2 

44.7 

-89.6 

40 

Test 

205.9 

144.7 

82.1 

56.5 

-53.0 


Posttest 

249.2 

143.7 

95.6 

40.4 

-98.5 


Pretest 

596.2 

124.0 

-286.5 

62.6 

-64.4 

12 

Test 

497.7 

82.6 

-294.7 

43.1 

-53.3 


Posttest 

549.1 

95.8 

-290.9 

62.9 

-71.2 


Pretest 

■ 


48.6 

-46.5 

124.0 

11 

Test 



79.9 

-98.4 

99.1 


Posttest 

152.6 

-108.3 

58.1 

-38.2 

131.8 


Pretest 

-530.9 

-220.4 

-152.2 

170.6 

-60.5 

36 

Test 

-452.6 

-266.1 

-91.7 

172.4 

-61.3 


Posttest 

-489.0 

-220.8 

-126.0 

134.6 

-63.9 


Pretest 


171.5 

118.9 

110.5 

-58.2 

37 

Test 


■ 

89.0 

-67.8 

-47.8 


Posttest 

-130.7 

1 mm 

101.3 

-105.3 

-60.9 


Pretest 

273.8 

-25.9 

-127.3 

-52.4 

-21.1 

18 

Test 

264.9 

-33.0 

-121.5 

-39.5 

-27.7 


Posttest 

309.1 

-30.7 

-126.5 

-48.9 

-17,7 


Pretest 

-299.0 

-17.9 

-104.4 

43.5 

-22.2 

8 

Test 

-292.1 

38.6 

-115.3 

25.7 

-23.9 


Posttest 

-302.4 

23.4 

-100.0 

70.5 | 

-23.9 
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TABLE 7 

Frequency Change and Damping of Solar Panel Modal Test 


Mode No. 

Pretest Math Model (Hz) 

Modal Test (Hz) 

Adjusted Posttest 
Math Model (Hz) 

Damping Ratio 

1 

37.0 

22.4 

22.4 

0.023 

2 

31.56 

26.8 

26.8 

0.017 

3 

34.55 

31.0 

31.0 

0.030 

4 

40.11 

- 


- 

5 

45.17 

38.5 

38.5 

0.030 

6 

59.70 

57.5 

59.70 

0.030 

7 

62.32 

64.9 

- 

- 


transformation, 



equation of motion, 



and force transformations, 

= t S ] Mt,A ' X ^T,A * 22e ^ 

i p| i s MWiM T> A 1 x} t,a (22f) 


^EE 

m t,a 


Wt,A t M ^T,A Mt,A 

(22h) 


B. Correlation 


The objective of the correlation of modal 
test data and analysis is to verify the mathe- 
matical model as shown in Eq. (22). However, 
this mathematical model is in a hybrid system, 
i.e., some degrees of freedom are expressed 
in the physical coordinates and others are in 
the generalized coordinates. The hybrid sys- 
tem is a result of the modal coupling technique 
used in the analysis. For a direct comparison, 
the solutions of Eq. (22) must be expressed in 
the form compatible to the measured test 
results. For this purpose, the normal modes 
obtained from Eq. (22) are transformed into 
the accelerometer locations. 

J T,A = C T 3 M t,A ^ 

where 

Mr . = analytical mode expressed 
’ at accelerometer locations 


The above values are the analytical predictions 
of the modal test configuration. 



normal mode solution of 
Eq. (22) 


The [m]t A term of Eq. (22d) can be 
expressed similarly to Eq. (21) as 


[M] 


T,A 


m t,a 

r:ER 

m t,a 


^re 

m t,a 

[^m ee ' 


T,AJ 


(22g) 


[T] = transformation matrix 

In Eq. (22), 1773 degrees of freedom exist in 
the hybrid coordinates, whereas [0 ']t A 
involves only 153 degrees of freedom. ’ Because 
the accelerometer locations often do not coin- 
cide with the analytical node points, Eq. (23) is 
an approximation. 
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For the modal test, a test mass matrix 
[m]t was constructed based upon physical 
mass distribution and corresponded to the 
experimental accelerometer measurements. 
Throughout the correlation, the analytical pre- 
dictions reduced to the accelerometer degrees 
of freedom [Eq. (23)] and the [m]t are used. 
The use of the [m] t ; A of Eq. (22) is compli- 
cated because it contains masses corresponding 
to hybrid coordinates and the degrees of free- 
dom of the test modes must be matched to cor- 
respond to the analytical degrees of freedom. 
The validity of [m] t is verified with the solu- 
tion of Eq. (22) by the mixed orthogonality 
check. 

Mt,a C m 1t E 0, Jt,a = [ m ']t,a 

Ideally [M r ]x a is a. diagonal matrix, however, 
as shown in Table 8, there are off-diagonal 
terms. In general for the first twelve modes, 
the off-diagonal terms are very small (less 
than 5%). This indicates that [M]t is indeed 
a valid mass matrix representing the total 
structural system for the first 12 modes. Addi- 
tionally the magnitude of the off-diagonal terms 
indicates the best accuracy one can expect in 
an orthogonality check of the test modes with 
the [M]t. Since a more detailed mass distri- 
bution was used in the analysis than in the 


modal test data reduction program, the orthog- 
onality of [M’] T a deteriorates for the higher 
frequency modes. 

The test modal data together with analyti- 
cal prediction [0']x A antl [M]t are read into 
the computer file for processing. The correla- 
tion work is automated. 

C. Mode Identification 

Prior to any correlation, the identification 
of test modes to the corresponding analytical 
modes is required. This task is achieved by 
the following criteria: 

(25 > 

t h 

Here the j test mode is checked with all 
the analytical modes [0’]x A* The test 
mode corresponds to the analytical mode 
related to the largest k th term in {m|t- A per- 
fect correlation exists when the k^ term is 
unity with all the other terms zero. After the 
identification of a test mode to an analytical 
mode, the correlation program calculates the 
necessary information required for a detailed 
comparison between the test and analysis. 

Table 9 shows the cross-orthogonality check 
between the test mode 701 and the first 30 ana- 
lytical modes. For this case, the test mode 701 
corresponds to the third analytical mode. 


TABLE 8 

Mixed Orthogonality — Model VIII 


4.35 4.40 

7.48 

7.83 

10.92 

13.36 

14.64 

17.95 

18.81 

23.42 

24.28 

OC 1 (J -m 

Frequency 

ZD . lo 

(Hz) 

1 2 

3 

4 

5 

6 

7 

8 

9 


11 



1 z -■ 


100 0.9 

-0.3 

-0.2 

0.8 

-0.9 

1.4 

1.6 

-1.7 

-0.3 

-0.9 

0.4 

l 

100 

-0.2 

-0.2 

0.3 

0.1 

-0.1 

0.2 

-3.1 

0,4 

-2.2 

-0.4 

2 


100 

1.3 

-1.6 

1.1 

-1.7 

-1.4 

0 

-2.1 

-1.2 

0.7 

3 



100 

-0.4 

1.2 

-0.2 

-0.3 

-2.7 

1.7 

-.03 

1.1 

4 




100 

1.0 

0.8 

1.1 

-0.6 

2.2 

-0.2 

1.1 

5 





100 

0.4 

-0.7 

0.8 

0,9 

-0.9 

1.8 

6 






100 

-1.8 

0 

0.2 

-2.3 

0.9 

7 







100 

-0.2 

-1.4 

1.2 

-0.2 

8 








100 

-0.9 

0.7 

-1.4 

9 









100 

-1.0 

0.2 

10 










100 

-2.3 

11 











100 

12 
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TABLE 9 

Cross Orthogonality - Orthogonality of Test 
Mode No. 701; Run Name DTA701 at 
Frequency 7.84 Hz with Respect to 
All Analytical Modes 


Analytical - 
Mode No. 

Freq 

Ortho 

■ 

l 

9.35 

- • oO 1 

2 

9 • 9q 

- *00^ 

3 

7 • Ma 

-.997 

9 

7.83 

- *0 1 ^ 

s 

10*92 

-•0Q3 

6 

13.36 

0 1 & 

7 

19.69 

. 0i)6 

8 

17.9& 

• 0 2 3 

9 

18 . 81 

• OOQ 

10 

23.^2 

•006 

1 l 

2 9 * 2a 

|01> 

12 

26*18 

- * 0 0 8 

13 

28.72 

• CQl 

19 

29.9a 

. QQ 1 

IS 

31*36 

- . 02Q 

1 6 

33. &9 

- . oQ9 

17 

3H.6fl 

. 009 

16 

3 S . 6 0 

-•000 

19 

36.95 

.01 1 

20 

38.93 

.002 

2 1 

39 . 1 1 

- #088 

22 

90 .s« 

-.0*1 

23 

92. OS 

- • 0 1 5 

29 

93.1S 

-• 0^3 

2 S 

95.32 

-•009 

26 

9 5 • 8 0 

« . 0 1 l 

27 

S 1 • «0 

.003 

26 

52. 9Q 

-.007 

29 

S3 . 1 5 

• OO 8 

30 

59. 9** 

- .0 1 0 

rokREL a T i u n Choi ct l 

ANALYTICAL MODE 3 r R EaUENcY /.MS 


D. Modal Test Equation 

The mathematical equation governing the 
modal test structural system is: 

jf| T e iwt = [M] T !uf T + [C] T li( T 

+ [K] T |u| t (26) 

where 

juL = displacement vector for each 
™ DOF associated with 

accelerometer measurement 

[m] t = test mass matrix 


[c] T = test damping matrix 

Nr = stiffness matrix 

| f = force vector from the shaker 

In contrast to the modal analysis as shown 
by the homogeneous Eq. (22), Eq. (26) is a 
forced response equation. The external har- 
monic excitation is provided by the shakers 
used in the modal test. The solution of Eq. (26) 
can be expressed in terms of generalized co- 
ordinates as follows: 

Wt-MtWt (27) 

where [$]t is the normal mode matrix, which 
is the eigenvectors of the test configuration 
measured at accelerometer locations. Substi- 
tution of Eq. (27) into Eq. (26) and pre multipli- 
cation by[$]£ result in 

[*]J )f| T e iwt = [-M J T lq| T + ["C J T |qf T 
+ ['Kj T !'ltT < 28 > 

where 

= E^Jt Wt 

generalized mass (28a) 

["Kj T = [0 ]t [K] t M x 

generalized stiffness (28b) 

Jt - C^]x W T 

generalized damping (28c) 

<% = v/K n /M„ 

n*k eigenvalue (28d) 

The accelerometer measurements obtained 
during the test are the response from the 
shaker force instead of the normal modes. The 
response is expressed as: 

!u| t = 0] T [*]* )n T e 1<Jt (29) 
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where 


u 2 /K 


+ i 2 P n w w„ 


When the shaker frequency lj is tuned to 
one of the natural frequencies u n , only one 
term in the [''H n (w) J dominates if p n , the per- 
centage of critical damping, is very small and 
all the natural frequencies are well separated. 
Thus the following approximation may be 
obtained: 


A relatively "clean" normal mode can be 
measured when damping is small and the nat- 
ural frequencies of the structure are well 
separated. In general, the damping of the 
structure is indeed very small but not all the 
natural frequencies are well separated. 

Hence, errors are expected when two normal 
mode frequencies are not well separated. 

E. Modal Test Output 

The output of the modal test is 

eigenvector, 

0] x (30a) 

eigenvalue, 

C~U 2 -J T (30b) 


of Eq. (28a). Equations (30a) through (30e) are 
compared to Eqs. (22a), (22b), and (22e) through 
(22h). 

F. Frequency 

The first data correlated are the natural 
frequencies. Table 10 lists the natural frequen- 
cies from analysis and corresponding test fre- 
quencies together with the mode description for 
the first twelve modes. Except for the eleventh 
mode, which is a scan platform mode, all the 
frequencies match fairly well. 

G. Orthogonality 

The generalized mass, Eq. (28a), is normal- 
ized to a unit matrix by proper normalization of 
each individual mode. Ideally, 

Wt-M? [M] t W t -M (31) 

Any errors in the mode shape measurement or 
mass data produce finite off-diagonal terms in 
the generalized mass matrix. If all off-diagonal 
terms are small, the measured normal modes 
are orthogonal to each other with respect to the 
mass matrix. Therefore, the orthogonality 
check serves as an indication of the accuracy 
of the measured test modes and test mass 
matrix. The orthogonality check of the first 
twelve test modes is in Table 11. The off- 
diagonal terms are indeed small and are within 
the 10% goal established for the test. The 
largest term occurs between the first and sec- 
ond modes. Since Table 10 indicates the natural 
frequencies of the first two modes to be almost 
identical, accurate mode shape measurements 
for these two modes are difficult to obtain. 


and modal force coefficients 

]P( T (30c) 

Evaluating the generalized mass matrix 
similarly to Eq. (22g), including the masses 
corresponding to the rigid body displacements, 
results in 



where 


H. Effective Mass (Ref. 7) 

In principle, the number of independent 
normal modes in a structural system is equal 
to the number of degrees of freedom. Obviously 
some of these modes are highly localized and 
of minor importance as far as the load analysis 
is concerned. With a limited number of modes 
obtained in the modal test, criteria are required 
to establish that the measured modes include 
the significant structural modes. For this pur- 
pose, the generalized rigid-body mass is used 
in the following way. In the analysis, the gen- 
eralized rigid -body mass is defined, as in 
Eq. (22g), as [ MRR] t a or 


[ M t\ = [" M J T ( 30e ) 
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TABLE 10 

Modal Frequency Comparison 


Mode 

Frequency (Hz) 

Description 

Analysis 

Test 

Error (%) 

1 

4.35 

4.51 

3.5 

Bending in X 

2 

4.40 

4.63 

5.0 

Bending in Y 

3 

7.48 

7.87 

5.0 

Lander in 9 ^ 

4 

7.83 

8.30 

5.7 

Lander in Y 

5 

10.92 

11.51 

5.1 

Lander in X 

6 

13.36 

14.09 

5.2 

Lander in 9 y 

7 

14.64 

15.35 

4.6 

Oxidizer tank in Z 

8 

17.95 

19.49 

7.9 

Fuel tank in Z 

9 

18.81 

19.83 

5.1 

Lander in 9 ^ 

10 

23.42 

24.85 

5.8 

Fuel tank in Y 

11 

26.18 

29.54 

11.4 

Scan platform in 0y 

12 

24.28 



26.49 

8.3 

Lander in Z 


TABLE 11 
Orthogonality 




3 

4 

5 

6 

7 

8 

9 

10 

11 

12 


BhI 

1 

100.0 6.2 

-0.2 

1.1 

-0.3 

1.1 

-2.3 

-0.6 

-1.7 

0.6 

-2.4 

0.0 


100,0 

0.1 

-1.2 

-4.1 

-3.0 

-0.9 

-2.5 

1.0 

• 

to 

3.4 

-1.5 



100.0 

0.4 

0.8 

1.6 

-0.2 

3.5 

-1.5 

-0.5 

-0.7 

0.4 




100.0 

1.0 

1.3 

1.2 

-0.1 

-1.8 

1.2 

1.1 

-0.5 

5 




100.0 

0.6 

0.8 

4.6 

-1.0 

0.2 

1.7 

-0.4 

6 





100.0 

0.4 

-0.5 

1.3 

4.4 

-0.6 

1.6 

7 






100.0 

-0.2 

-0.1 

3.6 

-0.1 

-1.3 

8 







100.0 

-1.1 

-1.9 

-0.2 

-1.5 

9 








100.0 

5.9 

-2.9 

2.7 

10 









100.0 

1.0 

-3.4 

11 










100.0 

2.5 

12 











100.0 
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From the modal test results of Eq. (3 Od) , 


(36) 


T 

[m^Jt = [Mr] [m] t [|u!r] (32b) 

An experimental generalized rigid-body-mass 
derived from [M^Ejrp is defined as 

where 

[m RE ] t = [lu| R ] T [M] t C*3t (34) 

[m ER ]t ■ [m RE £ (35) 

If all the experimental modes are obtained, 
it can be proved that 


[m^It = [m rr ] t * [m rr ] t a 


Another check is [m RR ] t must be equal to 
the rigid body inertia property of the test 
configuration. 


In the modal test, a limited number of 
measured modes was used to calculate the 
[mRR] t , now defined as effective mass. The 
comparison of [mRR] t and [m^R] t indicates 
whether thejnajor modes were effectively sur- 
veyed. If [MRR] t is close to [MRR] T , the 
major important modes with respect to the 
restrained point were obtained. Table 12 
shows the summary of the effective mass of 
the first twelve test modes and corresponding 
analytical modes. The first twelve modes rep- 
resent over 90% of the effective mass with 
respect to the base of the VTA/CTA. 


TABLE 12 

Effective Mass in Percentage 


Mode 

Mass 

X (%) 

Y (%) 

z (%) 

« x <%) 

(%) 

«z <%) 

Analysis 

Test 

96.42 

1.76 

0.01 

1.34 

85.20 

0.14 

89.51 

7.97 

0.03 

6.12 

78.23 

0.88 

Analysis 

1.67 

87.47 

0.0 

63.66 

1.49 


Test 

3.30 

86.49 

0.03 

64.45 

3.15 

mm 

Analysis 

0.95 

0.28 

0.0 

1.75 

0.02 

55.21 

Test 

0.89 

0.35 

0.01 

1.88 

0.02 

56.91 

Analysis 

0.06 

5.60 

0.03 

28.95 

0.01 

1.94 

Test 


4.81 

0.10 

27.51 

0.0 

1.66 

Analysis 

0.05 

0.06 

0.56 

0.10 

r 6.40 

19.58 

Test 

0.0 

0.11 

0.58 

0.19 

7.80 

20.43 

Analysis 

0.20 

0.06 

5.63 

■a 

4.77 

9.83 

Test 

0.38 

0.01 

6.67 


5.65 

7.39 

Analysis 

0.0 

0.01 

49.02 

0.01 

0.40 


Test 

0.01 

0.0 

51.80 

0.0 

0.63 


Analysis 

Test 

0.02 

0.01 

12.98 


0.04 

Kiiffl 

0.02 

0.0 

12.52 


0.0 


Analysis 

0.0 

0.06 

0.45 


0.0 


Test 

0.0 

0.05 

— . 

0.90 

0.02 

0.0 

0.09 
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TABLE 12 

Effective Mass in Percentage (contd) 



Mass 

Mode 

X(%) 

Y (%) 

z(%) 

* x <%) 

v*> 

Em 

Analysis 

10 

Test 

0.0 


3.15 

0.03 

0.0 

0.0 

0.0 


0.58 

0.03 

0.01 

0.04 

Analysis 

11 

Test 

0.0 

0.02 

13.60 

0.07 

0.0 

0.23 

0.0 

0.01 

11.57 

0.02 

0.0 

0.06 

Analysis 

12 

Test 

0.0 

m 

5.02 

0.06 

0.0 

0.0 

0.0 


22.17 

0.20 

0.0 

0.06 

Analysis 

99.37 

95.37 

90.45 

96.12 

98.33 

88.75 

Total 

Test 

94.39 

100.24 



95.96 



100.20 

95.89 

89.15 


I. Mode Shape 

The eigenvectors or mode shapes are 
important in the load analysis. The mode shape 
is expressed in the form of modal deflections 
of each degree of freedom (DOF). In the mathe- 
matical model, each DOF is assigned to a num- 
ber shown in Fig. 2, the ODTM modal test 
configuration. The direction of the modal 


SUPERSCRIPT C DENOTES CAPSULE LOCAL COORDINATE 
SUBSCRIPT S/C DENOTES SPACECRAFT COORDINATE SYSTEM 
NUMBER IN PARENTHESES IS A NODE NUMBER 


001 ) 


□ 


□ 


78.31 In. 


Z S/C 


deflection is indicated by the last digit in the 
DOF. Translation in X, Y, and Z are 1, 2, and 
3 respectively. For example, 1013 represents 
the motion of node 101 in the Z direction. The 
modal deflection of each DOF is expressed in 
their own local coordinate system, which is 
shown in Figs. 6 to 10. 

The maximum amplitudes of the modal 
deflection are normalized to unity. The analy- 
sis and test are compared and the difference at 
each DOF is expressed as a weighted deviation 
(WD) defined as 


WD = 


[O’] 


T,A 


-[*]■ 


normalized 


(37a) 


SUPERSCRIPT S DENOTES SCAN PLATFORM LOCAL COORDINATE 


S/C 


y S/C 


s/c 



Fig. 6 - Local coordinate of VLC 


Fig. 7 - Local coordinate of scan platform 
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SUPERSCRIPT B DENOTES BUS LOCAL COORDINATE 
NUMBERS IN PARENTHESES ARE NODE OR BAY NUMBERS 


SUPERSCRIPT P DENOTES PROPULSION LOCAL COORDINATE 
NUMBERS IN PARENTHESES ARE NODE NUMBERS 



Fig. 8 - Local coordinate of bus Fig. 10 - Local coordinate of propulsion module 


SUPERSCRIPT D DENOTES CABLE THROUGH LOCAL COORDINATE 
NUMBERS IN PARENTHESES ARE NODE NUMBERS 



Fig. 9 - Local coordinate of cable trough 


Also, for each mode the standard deviation 
(RSS) is calculated as 


RSS 



(37b) 


where 

N = number of DOF 
(WD). = weighted deviation of i** 1 DOF 

In Tables 13 and 14, a typical mode shape 
comparison of a mode and the summary com- 
parison of the first few modes, with only the 
important DOF and the standard deviations, 
are given. 


TABLE 13 

Modal Comparison — Analytical Mode 3 vs 
Experimental Mode 701 
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TABLE 13 

Modal Comparison — Analytical Mode 3 vs 
Experimental Mode 701 (contd) 


DOF 

Analysis 

— 

Test 

Weighted 

Deviation 

34 

-•3| 4-02 

- . 1 12-02 

-.202-02 

35 

♦ 968-02 

.570-02 

.399-02 

36 

♦7Q7-02 

.162-01 

-.9 1 3-02 

4 1 

-.914-02 

- . 189-01 

,972-0 2 

42 

. 3 1 4*00 

.299+00 

♦ 147-0 l 

S3 

.696-0 1 

.755-01 

-.593-02 

44 

-.895-03 

.263-02 

-.353-02 

45 

-. 1 15-01 

-.17 1 -02 

-.981-02 

46 

. 748-02 

.129-01 

-.538-02 

5 1 

-. 106*00 

-• 1 10*00 

.312-02 

52 

.353*00 

.325*00 

.281-0 l 

53 

. 140*00 

. 121+OC 

• 187-0 1 

ss 

-.923-03 

.282-02 

-.375-02 

55 

-. 143*01 

- .889-02 

-.544-02 

56 

.820-02 

. 133-01 

-.506*02 

6 l 

-.847-0 1 

- .966-0 l 

. 1 19-01 

62 

.403*00 

.322*00 

,804-01 

63 

* 135*00 

. 113*00 

.229-01 

64 

-.854-02 

-.347-02 

-.507-02 

65 

-.976-02 

-.824-02 

-.162-02 

66 

. 139-0 l 

.158-01 

188-02 

7 1 

-.309*00 

-.397+00 

.878-01 

mm 

*474*00 

.507*00 

-.328-0 l 

m m" 

*707-01 

. 137*00 

- .665-0 1 

m ■ ' 

-.364-02 

-.447-02 

.827-03 

K Wk 

.496-02 

-.124-01 

. 173-01 

mm: 

. 145-02 

. 146-01 

-.132-01 

e l 

-.56 1 *00 

- ,55 1*00 

-.980-02 

82 

• 399*00 

•366+00 

.327-01 

83 

.346-01 

•t. 368 — Q1 

.734-01 

04 

-.982-02 

-.216-02 

-.766-02 

85 

-. 1 19-01 

-.917-02 

-.270-0 2 

86 

. 1 1 5-0 1 

.155-0 1 

-.396-02 

9 1 

-.297*00 

-.291*00 

-.617-02 

92 

• 806 + 00 

.816*00 

-.967-02 

93 

-.756-0 l 

-.999-01 

.243-01 

94 

-.976-02 

-.669*02 

- • 307-02 

95 

-.121-02 

-.782-02 

•662-02 

96 

• 100-0) 

.520-02 

.482-02 

| 0 l 

« 166+00 

.167*00 

.195-0 l 

102 

.851 *00 

.931 *00 

-.800-01 

1 03 

-. 161*00 

-.8 19-01 

-.786-01 

104 

.816-03 

- ,629-q 2 

. 7 | 1-02 

105 

-.856-03 

.230-02 

-.315-02 

1 Q6 

-.208*02 

-.172*02 

-.351-03 

1 1 1 

-.949-01 

-• 106*00 

. 1 08-0 1 

| 12 

•864+00 

.824*00 

.4C’-Q1 

113 

-. 124*00 

-.312-02 

-. 12 1*00 

114 

-* 1 08-0 1 

-.352*03 

- . lO^-Q 1 

115 

. 195-01 

•227-01 

-.326-02 

1 1 * 

• 103-01 

•913*02 

. 1 l 8*02 

121 

-.409*00 

-.361*00 

- .408-0 1 

122 

.953*00 

• »cq*oi 

- .468*0 1 

1 23 

-.187*00 

-.985-01 

-.889-01 

124 

-. 166-02 

.239-02 

• • 4o 5* 02 

125 

-.475-02 

. 382-01 

-.429-0 1 

126 

.401-02 

.452-01 

*.4| 2-0 1 

| 3 1 

.287-01 

.327-01 

*.397-02 


TABLE 13 

Modal Comparison — Analytical Mode 3 vs 
Experimental Mode 701 (contd) 


DOF 

Analysis 

Test 

Weighted 

Deviation 

1 32 

. 100*01 

.951*00 

.491-01 

133 

- . 140*00 

-.11 1+00 

-.207-01 

134 

.323-02 

.5 J 3*02 

*. 191-02 

1 35 

120-01 

.656-02 

-.185-01 

1 36 

.915-02 

.225-01 

-. 133-0 1 

I 4 1 

.460*00 

.459*00 

.825-03 

142 

.790+00 

.691*00 

.994-01 

143 

-.268-01 

- .246-fl 1 

-.215-0 2 

144 


,662-02 

-.688-02 

145 

-. 106-01 

1 28-02 

-.935-02 

146 

.138-01 

. 195-01 

-.570-02 

1 5 1 

• 197*00 

.210+00 

129-01 

152 

.795+00 

■ 744*00 

.504-0 1 

153 

•500-0 1 

142*0 1 

.643-01 

154 

,308-02 

- . 193-03 

.327-02 

1 56 

.141*03 

-.9 1 1-02 

.926-02 

156 

.589-02 

. 158-Q 1 

-.994-02 

16 1 

-.101-01 

.228-01 

-.329*01 

162 

.704*00 

. 770*00 

. 142-0 1 

1 6 3 

.690-01 

• 2 15-01 

.476-01 

164 

-.731-03 

.296*02 

- . 369-02 

165 

-.757-02 

-.831-02 

.738-03 

166 

. 152-01 

♦ 1 1 1-01 

.408-02 

(OH 

- .464*00 

-.460+00 

-.392-02 

1012 

-.374*00 

- • 394*00 

.202-01 

1013 

151-01 

-.206-01 

.547-02 

1014 

•426-02 

.429*02 

-.389-04 

1015 

. 1 70-0 1 

.165-01 

,462-03 

1016 

.398-01 

.426-01 

-.279-02 

30 1 1 

♦470*00 

,467+OC 

.2j9-02 

3012 

.4 19*00 

.433*00 " 

- • 1 46-0 l 

3013 

.340-01 

.175-01 

. 1 6 5-Q 1 

3019 

.499-02 

.512-02 

126-03 

3015 

-.289-02 

-.343-02 

.54 1-03 

30 16 

. 164*01 

. 171*01 

-.742-03 

3031 

.458*00 

.452*00 

,608-02 

3032 

- 1 165*00 

-• 196*00 

.313-0 1 

3033 

-.191-01 

-.426-01 

.236-0 1 

3034 

.428-02 

.33 1-02 

.969-03 

3035 

-•267-02 

-.364-02 

,978-03 

3036 

• 160-01 

. 168-0 1 

-.798*03 

401 1 

.387*00 

* 3 7 S + OO 

. 121-01 

40 12 

-.704-0 1 

-.598-01 

- . 106-0 1 

401 3 

.756-02 

-.286-01 

.362-01 

4014 

.373-02 

.609*02 

-.31 7-02 

4015 

.289-02 

-.499*03 

.339-02 

4016 

. 181-01 

. 192*01 

-.11 2-02 

501 l 

• 55 1 +00 

.566+00 

148-01 

50 12 

.244*00 

• 269 + 00 

* .245-0 1 

5013 

* 1 1 5-0 l 

-.225*01 

.340-01 

5014 

•340-02 

•274-02 

.652*03 

5015 

-.223*02 

-.713*03 

-.151-02 

5016 

. 173-01 

.204*01 

- . 3 | 3-02 

201 1 

. 687 + QQ 

.616*00 

.701-01 

2012 

- . 186 + QO 

-.207*00 

•2| 7-01 

2013 

.204*00 

. 102+00 

• 102 + 00 

2014 

- . 425-02 

,293-0 2 

• • 7 18-02 

2015 

. 106-0 1 

.231-01 

126-01 
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TABLE 13 

Modal Comparison — Analytical Mode 3 vs 
Experimental Mode 701 (contd) 


DOF 

Analysis 

Test 

Weighted 

Deviation 

20 1 * 

* 14 7*01 

. 198-01 

-.516-02 

H02 1 

.673*00 

.743+00 

-.698-01 

4022 

* 120*00 

. 159*00 

- ,390-C 1 

*♦023 

-.756-QI 

-.324-01 

-.432-0 1 

*♦031 

.673*00 

.810*00 

-. 144*00 

4032 

.477-01 

.833-01 

-.356-0 1 

wma 

-.648-01 

-. 122+00 

.371-01 

4Q41 

.673*00 

.734+00 

-.603-0 1 

4042 

-.338-01 

- .465-0 l 

. 127-0 1 

4043 

- , 756-0 1 

-. 120 + 00 

.449-01 

60 1 1 

-.330*00 

-.266+00 

-.717-01 

60 1 2 

-.260*00 

- * 259*00 


60 1 3 

-.600-02 

.341-01 

. : ■ 


TABLE 13 

Modal Comparison — Analytical Mode 3 vs 
Experimental Mode 701 (contd) 


DOF 

Analysis 

Test 

' 

Weighted 

Deviation 

6021 

6022 

6023 

6031 

6032 

6033 
604 1 

6042 

6043 

.379*00 
-.848*00 
.600-01 
-. 190 + 00 
• 1 20*00 
.301-01 
.239*00 
.673*00 
- .379-01 

.390+00 
- .852+00 
-.395-01 
-• 135+00 
.120-01 
- .S96-0 1 
.267*00 
.649*00 
.125-01 

-.193-01 
.437-02 
.996-0 1 
-.551-01 
. 108 + 00 
•897-01 
- . 280-0 1 
.243-01 
-.503-01 

EXPERIMENTAL AT 60 
MAX ANALYTICAL AT 74 
RS“S ERROR* , .153 + 00 


TABLE 14 

Summary of Mode Shape Comparison 
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J. Local Kinetic Energy 

The generalized mass as shown in 
Eqs. (22h) and (28a) is expressed in the tenso- 
rial form as 

M if ** “ki (38) 

f =1 k=l 

If the normal modes are in the form of 
velocity, Mj?E is in the form of kinetic energy 

(see Eq. (5)). Then each term in Eq. (38) rep- 
resents the modal kinetic energy associated 
with that particular DOF. In essence, this is 
an itemized generalized mass. The comparison 
of each individual term in Eq. (38) between the 
test and analysis provides detailed information 
about the mass associated with each DOF and 
the modal amplitude of that DOF. This is 
especially valuable for those degrees of free- 
dom where the mass is questionable, such as 
the propellant tanks. The standard deviation 
(RSS) for the local kinetic energy comparison 
for each mode is calculated similarly to the 
mode shape comparison. In Table 15, a typical 
mode is selected for the kinetic energy com- 
parison. Also, in Table 16, a summary com- 
parison of the first few modes and the standard 
deviation of the kinetic energy comparison is 
given. 


TABLE 15 

Kinetic Energy Comparison - Analytical 
Mode 3 vs Test Mode 701 


DOF 

Analysis 

Experimental 

1 l 

*08 

•o’ 

12 

• 3 l 

♦ 2* 

13 

• oc 

.00 

1 4 

*co 

• 00 

IS 

♦ Co 

• q0 

16 

• 00 

• oo 

2 t 

• 20 

• 2 1 

22 

• 06 

• os 

23 

• 00 

• oo 

24 

• 00 

• oo 

25 

• 00 

• oo 

26 

■ [Uil I 

• 00 

31 

• 05 

♦ 05 

32 

■ 05 

• 0 3 

33 

• 00 

• oo 

34 

♦ GO 

• oo 

35 

• 00 

• oo 

36 

♦ 00 

• ol 

41 

• Go 

• 00 

mm 

• 06 

• 05 


TABLE 15 

Kinetic Energy Comparison - Analytical 
Mode 3 vs Test Mode 701 (contd) 


DOF 

Analysis 

Experimental 

43 

• oo 

• oo 

44 

• GO 

• 00 

45 

■ 00 

• 00 

46 

• Oo 

• oo 

51 

•Ol 

•0 i 

52 

• to 

• o 8 

53 

• 02 

• ol 

54 

*00 

• oo 

55 

•01 

• oo 

56 

• 00 

• 00 

61 

•01 

• 0 1 

62 

■ 1 1 

• 02 

63 

•01 

•ol 

64 

• 00 

• oo 

65 

• oo 

• rO 

66 

*00 

• c i 

7 1 

• 06 

• 0 8 

72 

• 13 

• 14 

73 

.00 

*0 1 

74 

« Oo 

• CiQ 

75 

• oo 

• ot 

76 

• Oo 

• 0 1 

81 

• 16 

• 14 

82 

• 06 

■ 06 

83 

• 00 

■ oo 

84 

• 00 

• GO 

85 

• 00 

• 00 

86 

• Oo 

• 01 

9 l 

• 09 

.08 

92 

• 66 

• 62 

93 

•01 

• 0 1 

94 

• Oo 

• oo 

95 

• 00 

• 00 

96 

•oo 

• oo 

101 

• 02 

• 0 2 

mm 

•50 

• 54 

■a 

• 02 

• QO 


• Oo 

■ oc 

105 

• 00 

• oo 

106 

• Oo 

• CO 

1 1 1 

•01 

• Cl 

112 

• 5 i 

• 42 

1 13 

•01 


1 14 

•01 

• 00 

115 

• Ol 

• 02 

116 

• Oo 

• 00 

121 

• 1 6 

• 1 1 

1 22 

• 84 

• 84 

123 

•03 

•Ol 

124 

• Oo 

• 00 

1 25 

• Oo 

• 0 6 

126 

• Oo 

♦ OS 

131 

• 00 

• oo 

132 

1 • 0 2 

• 84 

133 

• 02 

• Cl 

134 

• oo 

• oo 

135 

• Oo 

• oo 

1 36 

• oo 

• ol 

14 1 

• 12 

• 1 1 
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TABLE 15 

Kinetic Energy Comparison — Analytical 
Mode 3 vs Test Mode 701 (contd) 


DOF 

Analysis 

Experimental 

■n 

• 3 5 

• 2 4 


• Go 

• oo 


• 00 

• oo 


• Go 

• 00 


•00 

• o o 


• 02 

*0 2 


• 39 

• 3 1 

I S3 

• 00 

• oo 

154 

• 00 

• po 

1 55 

• 0Q 

• oo 

156 

*00 

• ol 

161 

• Q0 

• 00 

1 62 

•5 l 

• 45 

163 

«0c 

• oo 

164 

• Co 

• oo 

1 65 

• 0g 

* 0 0 

166 

•01 

• oo 

j 0 1 1 

5*39 

«.75 

10 12 

3.98 

4 * 03 

10 13 

•01 

•0* 

1014 

•43 

• 40 

1015 

7*59 

6*56 

1016 

62*2o 

6 4*90 

301 1 

4 * 5 & 

4 • J 0 

3012 

3 • 6 Q 

3 • 5 1 

3013 

♦ 03 

•0 1 

3014 

•02 

• 02 

3015 

♦Co 

• 00 

30 1 6 

• 1 1 

• 1 1 

303 1 

2.63 

2*34 

3032 

•33 

• 43 

3033 

•01 

• 0 2 

3034 

♦ 02 

• oi 

3035 

• 02 

• 0 3 

3036 

• 1 0 

• 10 

4011 

• 1 4 

•12 

40 l 2 

« 00 

« oo 

4013 

• 00 

• oo 

4014 

• 00 

• oo 


The local kinetic energy information is 
also valuable in identifying the mode shapes by 
describing the items with the largest kinetic 
energy. The data are also used in defining the 
system damping from the substructure damping 
measurements. 


TABLE 15 

Kinetic Energy Comparison - Analytical 
Mode 3 vs Test Mode 701 (contd) 


DOF 

Analysis 

Experimental 

4015 

• Oo 

• oo 

4016 

•03 

• 0 3 

501 1 

• 1 9 

• J 8 

5012 

• 04 

• 04 

5013 

• 00 

• 00 

5014 

• 00 

• 00 

5015 

• Co 

• 00 

5016 

• 02 

• 02 

2011 

• 9s 

• 70 

2012 

• 0 6 

• 10 

2013 

• 07 

• 0 3 

2014 

*•02 

• 0 2 

201 5 

•24 

• 60 

2016 

-•03 

• 09 

mmm 

• 05 

• 0 5 

1 

• oo 

• oo 

4023 

• oo 

• oo 

403 1 

• 05 

• 07 

4032 

• oo 

• oo 

4033 

• Oo 

•oo 

404 1 

. , ■ EfnSslali 

• 05 

4042 

• Oo 

• oo 

4043 

• oo 

• 00 

601 1 

• 02 

• 01 

60 1 2 

♦ 0| 

• oi 

6013 

♦ 00 

• oo 

6021 

• 02 

• 0 2 

6022 

• 1 0 

• 09 

6023 

• Oo 

• 00 

6031 

•01 

• 00 

6032 

• Oo 

• oo 

6033 

• Oo 

• oo 

604 1 

•01 

• 01 

6042 

• 07 

« 06 

■MB 

• Oq 

♦ pC 

R5S ERROR* *247400 


Since the magnitudes of the analytical modal 
forces are arbitrary, a normalization factor is 
found to multiply the analytical value to make a 
meaningful comparison. This normalization 
factor is calculated as follows: 


K. Modal Force 


In the analysis, the modal forces of the 
truss members are calculated for each mode. 
During the test, strain gauge data were taken 
at these truss members from which the test 
modal forces are obtained. The calculated 
modal forces from the analysis and measured 
member forces from the test are compared. 



(39) 
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TABLE 16 

Summary of Kinetic Energy Comparison 


^\Mode 

DOF 


1 

2 

3 

4 

5 

6 

Bus (%) 

A* 

11.45 

11.23 

6.87 

1.98 

13.76 

17.82 

T* 

11.38 

11.43 

6.25 

1.54 

11.90 

14.41 

Lander (%) 

A 

42.75 

19.22 

79.61 

77.85 

42.72 

43.65 

T 

40.58 

21.88 

80.65 

77.35 

” 

44.25 

42.51 

Oxidizer 


25.64 

39.35 

8.32 

11.16 

23.81 

21.05 

tank (%) 


26.42 

39.80 

7.75 

11.64 

22.51 

24.04 

Fuel 


15.26 

24.53 

3.10 

7.08 


7.66 

tank (%) 


16.01 

21.22 

2.93 

7.52 


9.18 

Pressure 

A 

0.84 

0.79 

0.18 

0.05 

0.23 

0.73 

tank (%) 


0.97 

0.78 

0.16 

0.05 

0.18 

0.86 

Engine (%) 


0.93 

2.00 

0.25 

1.02 

1.31 

0.84 

HB 

0.92 

1.78 

0.25 

1.14 

1.39 

0.85 

Scan 

A 

2.08 

1.92 

1.27 

0.64 

3.22 


platform (%) 

T 

2.01 

2.09 

1.63 

0.59 

3.80 

mmm 

RSS {%) 


0.452 

. J 

0.319 

0.247 

0.208 

0.233 

0.330 

♦Analysis. 
* Test. 









where 

a = normalization factor 
N = number of members 


together with the ratio of the comparison for a 
typical mode* The truss member identification 
is defined in Ref. (8). Only the axial forces are 
used in the comparison. The bending forces 
are not included in the comparison. 


P A = analytical modal force of the 
1 ith member 

T th 

p. = test modal force of the i member 


The factor a minimizes the difference 
between the analytical value and test value in 
the least-square sense. The analytical member 
forces are multiplied by a and compared with 
the test value. Table 17 shows this comparison 


The discrepancies between the test and 
analysis are thought to be from two sources. 
One is the inherent disadvantage of using strain 
gauges that were setup for high strain reading 
from the static test. Since low levels of strain 
were recorded in the modal test, the accuracy 
of the reading is in question. The other source 
of discrepancy is the difference between the 
analysis and test mode shape. The modal 
forces are very sensitive to local differences 
in the mode shape. 
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Member 

P (test) 
Mode 701 
Freq = 7.840 Hz 

r EC FRA Ml 


BUS658 

-•2319+02 

BUS 6 6 0 

-4773+01 

BUS662 

•2573+02 

BUS6S4 

- .9597+00 

l 

SPACECRAFT TRUSS 

BUSS36 

-.4991+02 

euS687 

.1795+02 

9 US 6 33 

.6392+02 

6US609 

- -11 10+03 

BUS63D 

-.5990+02 

EUS69 1 

.2747+01 

3US692 

-.1106+02 

BUS 69 3 

.6297+02 

B US634 

.2033+03 

BUS69 r 

-.1613+03 

3US695 

-.1849+02 

BUS697 

.6590+02 

1 

UPPER FLA f^E TRUSS 

BUST26 

-.5537+02 

BUS727 

.6018+02 

3 US 729 

-.9699+01 

B U S 7 3 C 

- .2076 + 02 

3 US 732 

.7323+01 

BUS 74 2 

.1023+03 

3 US 746 

-.1035+03 

lander truss 


BUS75C 

.1195+01 

3US751 

-.1671+03 

BUS 75 2 

.9908+02 

BUS753 

-.1791+03 

BUS754 

.2112+03 

3 U S 7 5 5 

.4500+02 

1 

BUS MAIN LONGERON 

BUS306 

.3042+02 

BUS910 

.1632*02 

BUS311 

-5996+01 

B US 31 3 

-.1721+02 

BUS 3 1 6 

- .1381+03 

B US 31 3 

-.1525+03 

BUS820 

24821+01 

BUS 323 

-.1869+02 

BUS 826 

.20 78 + 03 

9 US 930 

.1413+00 


TABLE 17 

Axial Force Comparison — Factor = 6.076 


P (analysis) 
Mode 3 

Freq = 7.480 Hz 


Ratio 

P (analysis) 


•- 27 36 + 02 
•5332+01 
. 3009+02 
•4137+tl 


-.4710+02 
. 1 738+02 
•6239+02 
■ • 1 C 96 + C 3 
-. 6524 + 02 
.2310+01 
.9632+01 
•6680+02 
. 1998+03 
•1533+02 
• 19 39+02 
.5352+02 


-.4972+02 
.6 221 + 02 
-.1255+02 
-.2627+02 
• 99S7+G1 


- . 92 94 + 02 


.6561+01 

■•1615+03 

•9439+02 

.1794+03 

7 

. 50 33+02 


.3631+02 
.1527+0? 
.1418+02 
■• 1518+02 


■•1313 + 03 
■ 5 1 56 + C 1 
.1806+02 
•2228+03 
.9394+01 



♦llBu+Cl 

. 1222+01 

.1167+01 

•4311+01 


. 9437+00 

. 1002+01 

• 9052 + 00 
•9787+00 

• 10 3 9+ Cl 
.3409+00 
. 3709+C0 


• 9305+00 
•9504+00 
•1C49+01 


. 9793+C0 
.1050+01 
.1294+01 

• 1 2 6 5 +G 1 
•1359+01 
•8096+00 

• 8557+00 


• 9667 + 00 
•9624+00 
. 1001+01 
•1055+01 


.1194+01 
. 9079+00 
.2365+01 
.3323+00 
.1097+01 
. 9611+00 
.1069+01 
. 9663+00 
.1072+01 
.5649+02 
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TABLE 17 

Axial Force Comparison — Factor = 6.076 (contd) 


Member 

P (test) 
Mode 701 
Freq = 7.840 Hz 

P (analysis) 
Mode 3 

Freq = 7.480 Hz 

Ratio 

P (analysis) 
P (test) 

BUS931 

- *3867+02 

-.3952+02 

.1022+01 

B US 93? 

- 5098+02 

.4794+0? 

. 2394+00 

EUSS35 

-.1038+03 

- .1059 + 03 

.1020+0 1 

8 US 339 

-.4313+01 

-.1583+01 

. 3625+00 

BUS 99 0 

.1417+02 

,1782+02 

.1257+01 

3 US 341 

-.2594+02 

--.2123+02 

. 9194+00 

PROPULSION S 

J3SYSTCV 



PPPQC 3 

.7637+01 

,3336+01 

,4369+ OL 

PPP004 

-.5652+01 

. 91 20+ Cl 

-.1614+01 

PPPOC 8 

.8962+01 

. 1 566+02 

.1747+01 

PPP011 

.1276+02 

.1552+02 

. 1216+01 

PFP012 

-.6929+02 

-.6563+02 

.96 1L+C0 

PPP019 

1603+01 

-.1349+02 

. 9391+01 

PPP028 

.9454+00 

.2602+01 

.2 753 +01 

PPP036 

.5470+02 

.4360*02 

.7970+00 

PFP037 

.2029+02 

.8557+01 

.4217+00 

PPP04C 

-.2429+02 

-.2961+02 

.1213+01 

FPPOM 1 

.1826+02 

.2142+02 

.1172+01 

P PP 043 

. 5946*01 

.6036+01 

.1015+01 

PFPC5? 

-.5330+01 

-.3619+01 

.6790+00 

P PPG92 

.2706+00 

-.2533+01 

-.9361+01 

R . S.S- 

E^ROR OF TEST VS ANALYS 

IS COMPARISON OF AXIAL F 

CRCES 



NORMALIZATION 

RSS FRRCR 



P ACTOR 

OF FORCES 

INPUT 

: ACT OR 

. 000 

.etc 

CALCULI 

k T ED FACTOR 

6.076 

6.937 

C .E *- 

.5462+00 



W.E .r 

.1632*00 




Several methods are used to establish a measure of error. They are: 


Description 


Relation 


Least square error (LSE) 



Reason 


No weight 
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Description 


Relation 


Reason 


Weighted least square 
error (WLSE) 



Weighted readings 
with higher 
strain 


Composite error (CE) 



Weighted readings 
with higher 
forces 


where 

= test strain gauge reading 
e „ = full scale strain gauge reading 


The results of the various weighting 
methods are in Table 18. Also the composite 
error is applied only to the Viking lander cap- 
sule adapter (VLCA) and V-S/C-A as shown in 
Table 19 to compare the data with the proof 
launch spacecraft Modal test (Ref. 9). 


TABLE 18 

Modal Force Comparison Errors 


Mode 

Least Square 
Error (lb) 

Weighted 
Least Square 
Error (lb) 

Composite 
Error (%) 

Analysis 

Test 

1 

708 

53.01 

2.09 

74.82 

2 

703 

32.28 

0.90 

50.62 

3 

701 

6.94 

0.94 

54.62 

4 

702 

17.69 

0.030 

19.40 

5 

704 

9.19 

0.015 

19.89 

6 

705 

6.38 

0.086 

21.94 

7 

711 

28.41 

1.16 

35.54 

8 

717 

42.66 

1.37 

75.13 

9 

712 

45.87 

1.73 

72.12 

10 

707 

2.26 

0.0035 

51.87 

11 

714 

22.39 

0.29 

101.3 

12 

713 

52.71 

1.61 

98.83 
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TABLE 19 

Composite Error of VLCA and V-S/C-A 


Mode 

VLCA 

(%) 

V-S/C-A 

(%> 

Analysis 

Test 

1 

708 

55.75 

46.81 

2 

703 

54.56 

41.67 

3 

701 

18.81 

12.06 

4 

702 

5.37 

7.19 

5 

704 

8.20 

17.52 

6 

705 

25.46 

8.53 

7 

711 

21.82 

14.89 

8 

717 

49.18 

51.88 

9 

712 

79.37 

82.52 

10 

707 

29.95 

59.40 

11 

714 

102.0 

108.6 

12 

713 

62.99 

21.99 


To help resolve the source of discrepancy, 
test force coefficients are checked using the 
modal acceleration data. The truss loads 
resulting from the modal acceleration times 
the mass applied as a static load are compared 
to the modal force coefficient. The compari- 
son for the VLCA members is shown in 
Table 20. The good correlation indicates the 
error to be the result of the difference in the 
analytical and test modes. 


TABLE 20 

Comparison of Modal Force by Inertia Load — 
Mode 713, Frequency 26.49 Hz 


Member 

No. 

Test 

(lb) 

Inertia 

db) 

Analysis 

db) 

750 

161 

175 

134 

751 

53 

47 

-20 

752 

262 

299 

166 

753 

279 

291 

290 

754 

116 

131 

-44 

755 

104 

130 

157 

LSE 


20.92 

85.41 

CE 


12.76 

62.99 




1 


Future efforts will be to obtain the force 
coefficient throughout the structure by applica- 
tion of the modal acceleration times the mass 
as an external load. 

L. Strain Energy 

Strain energy (SE) of a mode is a signifi- 
cant parameter in the description of a mode 
shape, especially when modal forces are of 
interest. Equation (5) directly relates strain 
energy to the stiffness matrix. 

For a mode, strain energy values describe 
the member or sets of members with the largest 
forces and deformations whereas kinetic 
energy (KE) describes the masses with the 
largest inertial motion. The significant dif- 
ference is the SE of the system is directly 
related to relative displacements at the ends 
of the members whereas KE is directly 
related to the absolute displacement. 

The axial strain energy of the i^ member 
monitored during the modal test is 


(SE). 


1 F? L. 

1 l i 

2 A~E~ 

l l 


(40) 


where 

(SE)^ = strain energy of the i^ member 

th 

F. = force in the i member 

L^ = length of the i^ member 

A. = cross-sectional area of the 
i^h member 

E. = modulus of elasticity of the 
ith member 


The SE of each axial member for both 
analysis and test is listed in the order of its 
magnitude. Also the SE of each member group 
is listed to help identify the characteristics of 
the mode shape. Table 21 shows the strain 
energy of each member and the comparison 
between test and analysis of a typical mode 
both in its magnitude and percentage. Table 22 
shows the comparison of the same mode listed 
in the order of magnitude. Table 23 shows the 
percentage SE of each member group for the 
first twelve modes. 
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TABLE 21 

Strain Energy Comparison - Analytical Mode 3 vs Test Mode 701 


Member 

Test (in. -lb) 

Analysis (in. -lb) 

Test (%) 

Analysis (%) 

: ENFRAME 





9US658 

. 4664-02 

.5497-02 

.2716-02 

.3975-02 

E U S 6 P r 

• 2 7 7T- 0 3 

.4 14 5-03 

.1616-03 

.2472-03 

B US66? 

. 5656-02 

.7705-02 

* 3 2 9 4-C 2 

.4596-02 

EUS6C4 

.1115-04 

. 2079-03 

.6511-05 

• 12 39-0 7 

TOTAL" 

. 1061-01 

.1432-01 

.6173-02 

.9342-02 

-p acecraft 

TRUSS 




BUSG9P 

.1 43 r -Cl 

.1275-01 

.9355-02 

.7621-02 

BUS6°7 

. 2537-02 

.2646-02 

• 1 c 3 6-G 2 

.1573-02 

EUS68S 

.3074-01 

.2519-01 

.1790-01 

.1503-01 

a ussr^ 

. 7540-01 

.7319-01 

.4449-01 

.4365-01 

EUSG90 

.3551-01 

.4213-01 

.2063-01 

.2513-01 

BUS591 

. 4129-04 


.2404-04 

.1741-04 

BUS69Z 

.6694-03 


.3898-03 

.3028-03 

B US693 

. 3192-01 


.1959-01 

.2136-01 

EUS6°4 

.2507+00 


■ ncHciil 

.1438 +00 

9USS°5 

. 15 4 P +00 

.1397+00 

.3003-01 

.8330-01 

GUS65F 

. 336 r -02 

. 3705-02 

.1960-02 

.2210-02 

3US897 

.2347-01 

.1351-01 

.1367-01 

.1104-01 

TOTALr 

-6244+CC 

. 595 2 + 00 

,3636+CG 

.3550+00 

UPPER PLATE 

TRUSS 




3 US726 

. 999 Q -D2 

.7741-02 

.5 922-02 

.4617-02 

BUS 72 7 

.1133-01 

. 1251-01 

.6600-02 

.7459-02 

B US 725 

.2151-03 

. 3603-03 

.1253-03 

.2149-03 

B U S 7 3 C 

.1613-02 

. 2583-02 

.5393-03 

.1541-02 

B US 73 2 

. 1233-03 

.2235-03 

.7209-04 

.1363-03 

BUS742 

.8196-01 

.6791-01 

.4773-01 

.4045-01 

3 US74P 

. 1056+00 

. 7729-01 

.6147-01 

.4610-01 

TOTALS 

.2108+00 

. 1685+00 

.122B+CC 

.1005+00 

LANDER tpls 

S 




B US 750 

. 1093-04 

.3265-03 

.6309-05 

.1S49-03 

B U S 7 r 1 

.2261+00 


.1317+00 

.1260+00 

3 US 752 

.5159-01 


.3004-01 

.2950-01 

BUS 75 3 

.1608+00 


.9365-01 

.9619-01 

BUS 754 

.3513+00 

.4013+00 

.2104+00 

.2397+00 

B'JS7 cr 

.1647-01 

.2101-01 

.5592-02 

.1253-01 

TOTAL- 

.9153+00 . 

. 3435 +GC 

• 4 7 5 3 + 00 

.5031+00 

BUS MAIN LO 

NG r RCN 




BUSSCf 

*2185-03 

.3113-03 

. 1272-03 

.1857-03 

B US 310 

.4136-03 

.3409-03 

.2409-03 

•2G33-03 

BUSS1 1 

.1798-04 

.9997-04 

. 1041-04 

.5963-04 

B US 313 

. 5177-03 

.4030-03 

.3015-03 

.2404-03 

BUS316 

.4413-02 

.5312-02 

.2570-02 

.3169-02 

B US 31 3 

. 5333-02 

.3995-02 

.3133-02 

.2383-02 

BUS 92 0 

.3215-04 

. 3677-04 

.1972-04 

.2193-04 

B US 923 

.6257-03 

.5843-03 

.3644-03 

.3435-03 

BUS 92 F 

.1033-01 

. 1187-01 

.6C 15-02 

.7081-0? 
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TABLE 21 

Strain Energy Comparison - Analytical Mode 3 vs Test Mode 701 (contd) 


Member 


C us 3 30 
BUS 3 3 1 
3 US 33? 
BUS 33 c 
3 US 33 ^ 

Eussur 

3 US 341 

TOTAL r 


Test (in. -lb) 


. 3011-07 
.3941-02 
.4339-02 
.2546-02 
. 2564-04 
.4 8L2-C3 
.1170-02 

. 3 4 9 f - 0 1 


ROFULSICN SUBSYSTEM 


P °P0G3 
PPP004 
POPOO? 
FTPD1 1 
PPP01? 
PPFQ1? 
PPPQ29 
PFP03P 
P PP 037 
prp04t 
PPP041 
PPF04 7 
PPP053 
FPP08 2 

TOTAL- 

TCTaL= 


. 1049-03 
.6654-04 
.4594-03 
.3571-03 
.5433-02 
.3424-05 
.2150-05 
• 6 6 3 3 - O 1 2 
. 794G-C3 
.1057-02 
. 7524-03 
.2331-03 
.2297-03 
.2134-00 

. 2C19-C1 

.1717*01 


INPUT FACTOR 
CALCULATES FACTCP 


Analysis (in. -lb) 

Test {%) 

Analysis (%) 

.1331-03 

.1753-07 

.7937-04 

.4116-02 

•2295-02 

.2455-02 

.4262-02 

.2 913-G2 

.2542-02 

. 2650-02 

, 1 4 8 3 - C i 2 

.1591-0? 

.3500-05 

.1551-04 

.2033-0^ 

.7590-03 

.2796-03 

.4527-03 

.7335-03 

.6312-03 

.4673-03 

. 3566-01 

.2036-01 

,2127-Ll 

.2002-04 

.6109-04 

.1194-04 

.1733-03 

.3875-04 

.1034-03 

.1402-02 

.2675-03 

.3363-03 

.5290-03 

•2080-03 

,.3149-03 

.3711-02 

.5493-02 

.5196-0? 

.6931-03 

• 4906-D 5 

.3539-0? 

.1630-04 

.1252-05 

.9721-06 

.4213-02 

.3863-02 

.2513-02 

.1413-03 

.4627-03 

.3430-04 

. 1570-02 

.6154-03 

.9365-0? 

.1036-02 

.4392-03 

.6177-03 

.2917-03 

.1649-03 

.1740-03 

.1054-03 

.1 332-03 

.6239-04 

. 1B7C-L4 

.1243-06 

.1115-04 

.1332-01 

.1175-01 

.1123-01 

. 1677+01 

.1000+01 

.1000+01 

VS ANALYSIS CCNFAFISC*/ OF STRAIN 5 

NERCY 

NORMALIZATION 

RSS ERROR 


FACTOR 

OP FORCES 


6.076 

6.937 


.933 

.001 



TABLE 22 

Strain Energy in Order of Magnitude - Analytical Mode 3 vs Test Mode 701 


Member 

Strain Energy (Test) 
(in. -lb) 

Member 

Strain Energy (Analysis) 
(in. -lb) 

BUS? 54 

.2104+00 

3US754 

.2397+00 

BUS £94 

.1460+00 

BUS694 

. 1430 + CC 

9US751 

.1317+00 

9 US 7 51 

.1260+00 

E US 753 

.9365-01 

EUS753 

. 9619-01 

SUS695 

. 9003-01 

BUSS 95 

.9330-01 

BUS 74 6 

.6147-01 

BUS746 

►4610-01 

9U S742 

.4773-01 

8USGB9 

.4385-01 

5 US63B 

.4449-01 

8US742 

.4045-01 

BUS752 

. 3 004 - Cl 

BUS 7 5 2 

.2 350-01 

5 US 69 t 

.2069-01 

6US69C 

•2513-01 
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TABLE 22 

Strain Energy in Order of Magnitude - Analytical Mode 3 vs Test Mode 701 (contd) 



M. Reaction Forc es 

The calculation of reaction forces at the 
base of the VTA/CTA by the modal forces 
from the V-S/C-A strain gages and the modal 
accelerations provides a check on the test data. 

The [MRE] T matrix from Eq. (30d) repre- 
sents the reaction forces from the modal 
accelerations. The reaction forces calculated 
from the V-S/C-A jnodal member force data 
are compared to [mREJ^ The comparison of 
the data are in Table 24. 


N. Generalized Mass from Modal Damping 

The generalized mass terms [m E ?J from 
Eq. (30d) are checked by the following relation- 
ship. At a normal mode, the damping values 
are offset by the force values. From Eq. (28), 


rcO T |i} T = [*]T{f} T (41) 

for the n mode. Since 
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Mode 

Bedframe 

<%> 

- A* 

0.01 

T* 

0.01 

2 A 

0.09 

T 

0.06 

3 A 

0.09 

T 

0.06 

4 A 

0.10 

T 

0.07 

5 A 

0.0 

T 

0.0 

6 A 

0.05 

T 

0.04 

7 A 

0.26 

T 

0.20 

8 A 

0.36 

T 

0.26 

9 A 

0.20 

T 

0.20 

10 A 

0.01 

T 

0.04 

n A 

0.10 

T 

0.04 

i-i 

to 

> 

0.06 

T 

0.10 


* 

Analysis. 

*Test. 


TABLE 23 

Strain Energy Distribution Comparison 







v - s/c 

Truss 

(%) 

Upper 

Plane 

Truss 

<%) 

Lander 

Truss 

(%) 

Bus Main 
Longeron 
<%) 

Propulsion 

Subsystem 

(%) 

0.86 

0.03 

0.09 

0.04 

0.02 

0.84 

0.03 

0.08 

0.04 

0.03 

0.73 

0.03 

0.04 

0.01 

0.10 

0.78 

0.02 

0.04 

0.01 

0.09 

0.36 

0.10 

0.50 

0.02 

0.01 

0.36 

0.12 

0.48 

0.02 

0.01 

0.35 

0.03 

0.42 

0.04 

0.06 

0.40 

0.03 

0.40 

0.04 

0.06 

0.54 

0.03 

0.30 

0.04 

0.09 

0.55 

0.04 

0.26 

0.04 

0.12 

0.43 

0.05 

0.31 

0.04 

0.13 

0.47 

0.07 

0.24 

0.04 

0.15 

0.22 

0.05 

0.11 

0.06 

0.35 

0.24 

0.08 

0.11 

0.07 

0.38 

0.09 

0.01 

0.09 

0.02 

0.44 

0.11 

0.02 

0.25 

0.03 

0.34 

0.05 

0.05 

0.44 

0.03 

0.23 

0.05 

0.04 

0.38 

0.04 

0.29 

0.08 

0.02 

0.57 

0.02 

0.30 

0.07 

0.02 

0.58 

0.03 

0.28 

0.22 

0.06 

0.46 

0.02 

0.13 

0.22 

0.27 

0.58 

0.01 

0.08 

0.08 

0.10 

0.63 

0,02 

0.11 

0.23 

0.03 

0.51 

0.03 

0.12 
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TABLE 24 

Comparison of Reaction Force 


Mode 

Reaction Force: Accelerometer (Strain Gage) a 

F 

X 

(lb) 

HH 

F 

Z 

(lb) 

M x 

(in. -lb) 

M y 

(in. -lb) 


708 

401 


-8 

16700 

60000 

1 


(517) 


(54) 

(15587) 

(52441) 


703 

75 

383 

8 

-52400 

11600 

1394 


(78) 

(391) 

(-102) 

(-33728) 

(6802) 

(1926) 

701 

45 

-28 

-4 

10300 

1110 



(50) 

(-36) 

(-10) 

(9540) 

(-1398) 

(12847) 

702 

19 

131 

19 

-49300 

59 

2440 


(14) 

(138) 

(-22) 

(-44093) 

(-1311) 

(2374) 

704 

-2 

19 

-44 

-3900 

23700 

8060 


(-D 

(-1) 

(-50) 

(-1585) 

(25842) 

(9087) 

705 

25 

-4 

-107 

684 

15000 

-3440 


(22) 

(5) 

(-119) 

(62) 

(14108) 

(-4483) 

711 

12 

-2 

1080 

241 

18100 

-2240 


(5) 

(-2) 

(1106) 

(974) 

(18207) 

(-3388) 

717 

17 

-6 

444 

1810 

833 

-1490 


(28) 

(-40) 

(465) 

(3392) 

(762) 

(-557) 

712 

2 

28 

125 

2600 

-438 



(17) 

(81) 

(97) 

(-6375) 

(248) 

(920) 


Numbers in parentheses refer to strain gage measurements. 


Each generalized mass term M nn is calculated 
from Eq. (43) by measuring the modal damp- 
ing, frequency, velocity of the participation 
factor, shaker force, and displacement in the 
force direction. The comparison is in Table 25. 
Since the modal damping p n varied, the M n n 
varied accordingly. As noted, the M nn is 
within the range of the damping measured; how- 
ever, because of the variation in p n , M nn can- 
not be accurately measured by this method. 


then 


C = 2 p M co 
nn r n nn n 


(42) 


M 


L # «lJ 


nn 2p u> 6 
K n n H n 


w 


(43) 
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TABLE 25 

Generalized Mass Comparison 


Mode 

Based on 
Damping 

Based on 
Mode Shape 

— 

High 

Aver- 

age 

Low 

708 

0.065 

0.141 

1.021 

0.261 

703 

0.111 

0.179 

0.617 

0.252 

701 

0.130 

0.227 

0.833 

0.290 

702 

0.063 

0.091 

0.410 

0.364 

704 

0.133 

0.290 

0.785 

0.344 

705 

0.056 

0.081 

0.330 

0.244 

711 

0.902 

1.499 

3,037 

0.881 

712 

1.111 

2.043 

5.557 

0.776 

707 

0.011 

0.024 

0.060 

0.083 

714 

0.171 

0.376 

3.141 

0.397 

713 

1.292 

2.179 

6.692 

1.101 


FINAL DYNAMIC MODEL 


The final VO dynamic model is a combina- 
tion of the modal test configuration and sub- 
structure characteristics verified by substruc- 
ture tests. Since Eq. (22) representing the 
modal test configuration was verified by test, 
the parameters are accurate. Modifications 
are made to update the parameters to the 
V-S/C Mission B configuration. The subscript 
FB denotes Mission B configuration. The 
changes are: 


CHANGES REASON 


[ M ] X;A t0 
Mp,T t0 Mp,FB 
Ms,T t0 Ws,FB 


The mass of the pro- 
pellant is changed to 
the Mission B config- 
uration; see Table 2. 

The dynamic charac- 
teristics of the 
propulsion module 
changes 

The dynamic charac- 
teristic of the scan 
platform changes 
since its joint is 
allowed to move along 
the serrations 


The parameters verified by the modal test 
and updated to represent Mission B are com- 
bined with the remainder of the substructures. 
Each substructure is verified by a modal test. 
The combined equation of motion represented 
by Eq. (18) becomes 


Mr 


( 

Mr 


Mr 

Mb 


Mb 


Mb 

Mp 


Mp 


Mp 

Ms 


Ms 


Ms 

Me 


Me 

CO 

Me 

falsPA 


Hspa 


M'spa 

Mspi 


Mspi 


Mspi 

MsP2 


Msp2 


MsP2 

MsP3 

FB 

MsP3 

k > 

FB 

MsP3 


Mb ■ [S]{u} B (45) 

Ml-tSlWiMi (46) 

Table (26) summarizes the various param- 
eters verified by the corresponding tests. 
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TABLE 26 

Test Verified Models 


Parameter 

Description 

Test for Verification 

[ m ] fb 

Mission B mass matrix 

All the modal tests and the propellant 
effective weight tests 

M 

Rigid body mode 

pjpj 

Calculations M FB should check with the rigid 
inertia properties of the V-S/C/CTA/CTA at 
its base 

N 

Model of bus with CT A/VTA 

Static test on the bus and system modal 
test; the CTA/VTA only included in the 
system modal test 

{ q fp 

Displacement function of 
propulsion 

Modal and static test of propulsion module 
and system modal test with referee fluid 

Ws 

Displacement function of scan 
platform 

Modal Test of scan platform and system 
modal test without joint slippage; the joint 
slippage was included for load analysis 
since higher forces would result in joint 
slippage 

He 

Displacement function of 
cable trough 

Modal Test of cable trough and system 
modal test 

WsPA 

Displacement function of solar 
panel with relay antenna 

Modal and static solar panel tests with 
a relay antenna 

{ q }sPi 
(i = 1, 3) 

Displacement function of solar 
panel 

Based upon a mathematical model of the 
solar panel with a relay antenna; the relay 
antenna was removed from the math model 
for this configuration 

( 1 

0 

1 1 

w 

Damping matrix for Mission B 
configuration 

Data based upon modal tests where avail- 
able; the damping matrix was diagonalized 
at each transformation; the kinetic energy 
evaluation of the modes was used as a 
guide to estimate damping; solar panel 
viscous dampers were not included but 
estimated as a modal damping 

WpB 

Stiffness matrix for Mission B 
configuration 

Substructure modal and static tests; sys- 
tem modal tests 

[s] 

Load matrix 

Substructure modal and static tests; system 
modal tests 
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(49) 


The real eigenvalues and eigenvectors of 
Eq. (44) with [C]pg = 0, results in 



Substitution of Eq. (47) with Eq. (44) and pre- 
multiplication of [$]£_ results in 

r D 




{ P }b ' [*] FB { x j FB 

(50 > 

The procedure to obtain the flight model 
was not complicated. The effort required to 
obtain was based upon engineering 

judgement. Since damping for most substruc- 
tures is measured or estimated, the system 
damping matrix was based upon the kinetic 
energy contribution of the substructures to the 
system mode. 

The data represented by Eqs. (48) through 
(50) was transmitted to MMA to couple the VLC 
on the VO to obtain the V-S/C mathematical 
model. 

CONCLUSION 

Highlights of the steps leading to the VO 
dynamic model are summarized with the 
mathematical equations and data. The success 
is attributed to substructure modal coupling 
concepts where each substructure is experi- 
mentally verified. The degree of correlation 
is dependent on the use of the dynamic model. 
The mode shapes and frequencies of the sys- 
tem correlated well but some difficulty exists 
in correlating modal forces. A measure to 
establish the degree of correlation is still 
required. 

An overall structures and dynamics pro- 
gram integrating the substructure analysis, 
hardware availability, and test will result in 
a successful dynamic model. Overall system 
modal tests may not be required. 
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